Computational Screening of the Phytocompounds from the Plant Ballota nigra Linn against the Human Papillomavirus (HPV) E6
Nishandhini Marimuthu1, Viswanathan T1, Mahendran Radha*, Jeyabaskar Suganya
Department of Bioinformatics, School of Life Sciences, Vels University, Pallavaram, Chennai-600117,
Tamil Nadu, India.
*Corresponding Author E-mail: mahen.radha@gmail.com
ABSTRACT:
Cancer is major cause of death rate in world wide. Cervical cancer is the second most common type of cancer in women, where the cancer develops at lower part of the uterus. Literature studies clearly depicts that the mainstream of the human cervical cancers are caused by high risk HPV infections. Human Papillomavirus (HPV) E6 protein was the most prevalent high risk oncogene causing cervical cancer and still there is no proficient drug against this disease. Many researchers are involved to combat this cancer using phytodrugs. The current research was to identify the potential plant based drug against the cervical cancer using docking studies. Docking studies of eighteen phyto compounds from the plant Ballota nigra of the family Lamiaceae were carried out using Argus Lab. From the results, Ballotinone and Ballonigrin were found to have a best binding efficiency with the active site residues of the protein. The binding energies of the 2 compounds were -12.0691 Kcal/mol and -10.2564 Kcal/mol. This study provides promising anti-cervial cancer inhibitor for further drug development.
KEYWORDS: Ballota nigra Linn, Phytocompounds, Binding energies, docking studies and Cervical cancer.
INTRODUCTION:
Cancer of the uterine cervix is second most common cause of female malignancies in worldwide1. Human papilloma viruses (HPVs) of small DNA that belong to papilloma virus family, which infect mucosal epithethelial cells, genital tissues and upper respiratory tract2. About, 99 % of the cervical cancer is caused by human papillomavirus3. According to the progression of malignancy, genital HPV types are divided into low-risk and high risk groups. Low-risk HPV types 6 and11 cause genital warts, which do not progress towards cancer. High risk HPV types cause malignant lesions that will promote to invasive carcinomas. High risk HPV types are seen in over 95 % of the cervical cancer4. Persistent HPV of high risk is found in cervical intraepithelial carcinomas5.
The HPV-16 E6 and E7 oncoproteins have ability to inactivate p53 and members of the retinoblastoma (Rb) tumour suppressor family respectively which play a major role in tumour formation. High risk viruses have been exposed to alter pathways in cell cycle and signal transduction pathways6.
HPV 16 E6 is formed by approximately 151 amino acid protein with two zinc finger joined by an interdomain linker of 36 amino acid, flanked by short amino (N) and carboxyl(C) terminal domains7. Inactivation of Mdm2 pathway, occurs when E6 protein arrest the proapoptoptic activity of p53 by interacting with a plethora of cellular proteins8. Ubiquitin mediated degradation occurs when E6 binds to p53 in the cytoplasm by blocking its translocation to the nucleus. Thus by preventing the p53 translocation will repress the tumour growth9.
Effective treatment against the HPV infection is still unavailable. So there is an urgent need for the promising drug against cancer. For many decades, Plant has been an important source of herbal medicine, because of the presence of various health’s sustaining constituents. The foremost objective of this research study is to identify the bioactive compound from the natural plants that produce physiological actions on the human body.
Ballota nigra (Black horehound) is a perennial herb of the family Lamiaceae. The plant is used to treat nervous disorders, stomach cramps, whooping cough, diabetes and gastrointestinal disorders10-12. Therefore, the present study was made to screen the eighteen compounds from the plant and their interaction with the active site residues of the protein HPV16 E6 using docking studies.
MATERIALS AND METHODS:
Computational Methodology:
Protein Preparation and energy minimization:
The three dimensional structure for Human Papilloma virus E6 oncoprotein (1VZN) was retrieved from Protein Data Bank. Then the protein chosen for the study should be optimized by removing water molecules, ions, ligands from the structure file. Further, the hydrogen atoms present in the protein are optimized and minimized according to Hartree-Fock calculation method of Argus lab.
Preparation of Ligand structures:
The natural compounds used for the docking study were selected from the literature studies. Eighteen phyto compounds from the plant Ballota nigra were chosen for the study. These structures were retrieved from Pubchem and some of the ligand structures were drawn using chemsketch. Further, optimization of the compound structure was carried out using Hartree–Fock (HF) calculation method using ArgusLab.
Active site Prediction:
Identification of active site is often the initial point in protein–ligand binding site prediction for structure-based drug design. Active site of the protein was obtained using the metapocket. Metapocket uses interaction energy and a simple Vander Waals probe to locate energetically favorable binding sites of the protein. The most potential binding site of the protein was analysed.
Molecular Docking Study:
Molecular Docking of the protein and ligands were performed using Argus Lab 4.0.1. Docking was performed using “Argus Dock” with exhaustive search docking function of ArgusLab with a grid resolution of 0.40 Å. Protein structure of HPV 16 E6 (1vzn) was loaded and the binding site of the protein were analyzed. Then the compounds chosen for the study were introduced and the docking analysis were performed by shape based algorithm and the flexible docking of the ligands by Ascore scoring function. Docking of the protein- ligand complex was evaluated by scoring function. 10 independent runs were conducted and one pose for each complex structure was obtained. The stability of the docking complex was analyzed by lowest Ascore and hydrogen bond interactions.
RESULTS AND DISCUSSION:
Natural inhibitors from the plant were docked with active site residues of the protein HPV 16 E6. The active site residues obtained from metapocket were listed in the Table1.Among the eighteen compounds, Ballotinone and Ballonigrin showed a lowest binding energy of -12.0691 Kcal/mol and -10.2564 Kcal/mol. Ballotinone was found to interact with E6 residues of Leu 45 and Glu 120 by forming hydrogen bonds, Where as the compound Ballonigrin was found to interact with the residue of Tyr 83 by forming hydrogen bond are shown in the Figure 1.Of all compounds, five compounds showed least binding energy and their interactions with the residue, hydrogen bond distance are given in the Table 2.
Table 1: Catalytic site residues of HPV 16 E6 protein
|
Amino acid residue in the active site |
Structural unit of the residues |
|
Gln13,Cys 37,Val 38,Cys 40,Lys 41gln 42,Leu 44,Leu 45,Glu 48,Val 49, Asp 51,Phe 52,Arg 55,Ser 78,Lys 79,Ser 81,Glu 82,Tyr 83,His 85,Tyr 86,Ser 89, Arg 109,Cys 110, Pro 116,Leu 117,Cys 118,Pro 119,Glu 120,Lys 122 and Leu 126 |
Helix |
|
Leu 19, Pro 20, Gln 21, Leu 22,cys 23,Thr 24, Leu 26, Lys101Cys 113, Gln 114and Lys 115 |
Coil |
|
Gln 97, Gln 98, Tyr 99, Asn 100, Leu 103and Asp 105 |
sheet |
a. Ballotinone with HPV 16 E6
b. Ballonigrin with HPV 16 E6
Figure 1: Docking of HPV 16 E6 protein with Ballotinone and Ballonigrin
Table 2: Interactions of HPV 16 E6 protein with compounds from Ballota nigra
|
S. No |
Compounds |
No of H-bonds |
Amino acid residue atom |
Ligand atom |
H-bond distance (Å) |
Binding energy Kcal/mol |
|
1 |
Ballotinone |
2 |
Leu 45 (NH) Glu 120 (NH) |
(O2) (O2) |
2.75 3.00 |
-12.0691 |
|
2 |
Ballonigrin |
1 |
Tyr 83 (OH) |
(O2) |
2.36 |
-10.2564 |
|
3 |
Caffeoylmalic acid |
3 |
Gln 98 (NH) Tyr 99 (NH) Tyr 86 (OH) |
(O2) (O2) (O2) |
2.48 2.53 2.61 |
-9.09384 |
|
4. |
Chlorogenic acid |
3 |
Tyr 83 (OH) Asp 32 (OH) Cys 37 (S) |
(O2) (O2) (O2) |
2.90 2.90 2.78 |
-8.5403 |
|
5. |
martynoside |
1 |
Lys 41(OH) |
(O2) |
2.52 |
-7.9880 |
Our docking result revealed, the compounds from the plant were known to inhibit he protein HPV 16 E6 of cervical cancer. Interaction of the compounds with amino acid residues of p53 binding, which prevent E6 protein from binding with p53 protein.
CONCLUSION:
The molecular interaction between the HPV 16 E 6 protein and p53 show a major crucial role in cervical cancer. Hence, molecular docking studies of the compounds against the active site of the protein HPV 16 E6. From the results, it is revealed that Ballotinone and Ballonigrin showed high binding affinity towards the protein by disrupting the human tumor suppressor protein p53-HPV oncoprotein complex. In conclusion, few compounds of the plant Ballota nigra could be the successful drug against cervical cancer. This study is further investigated in drug designing for development of effective therapeutic medicine against cervical carcinoma.
REFERENCES:
1. zur Hausen H. Papillomaviruses and cancer: from basic studies to clinical application. Nat Rev Cancer. 2002; 2: 342–350.
2. zur Hausen H, de Villiers EM. Human papillomaviruses. Annu Rev Microbiol 1994; 48: 427–47.
3. J.M.M. Walboomers, M.V. Jacobs, M.M. Manos, F.X. Bosch, J.A. Kummer, K.V. Shah, et al. Journal of Pathology.1999; 189: 12
4. Trottier, H and Franco, EL. The epidemiology of genital human papillomavirus infection. Vaccine, 2006; 24: S1-15.
5. Wright TC, Kurman RJ. A critical review of the morphologic classification systems of preinvasive lesions of the cervix: the scientific basis for shifting the paradigm. Papillomavirus Rep 1994; 5:175–82.
6. Dell.G and Gaston.K. Contributions in the domain of cancer research: Review, Human papillomaviruses and their role in cervical cancer CMLS, Cell. Mol.Life Sci.2001.
7. Howie, HL., Katzenellenbogen, RA. and Galloway DA. Papillomavirus E6 proteins. Virology, 2009; 384: 2, 324-334.
8. Thomas, MC. and Chiang, CM. E6 oncoprotein represses p53-dependent gene activation via inhibition of protein acetylation independently of inducing p53 degradation. Molecular Cell, 2005; 17: 2: 251-264.
9. Clinton Jones Cervical Cancer: Is Herpes Simplex Virus Type II a Cofactor Clinical Microbiology Reviews 1995.
10. PDR for herbal medicines, Medical Economic Co. Montvale, New Jersey.2000; 98-99.
11. Gunther RT. The Greek Herbal of Dioscorides, Hafner Publishing Co, New York. 1959; 347.
12. Darbour N, Baltassat F and Raynaud J. Sur la presenced’un O-heteroside et d’un C-heterosided’apigenindans les feuilles de Ballota foetida Lam .(Labiees). Pharmazie.1986; 41:
Received on 08.03.2017 Modified on 22.05.2017
Accepted on 12.06.2017 © RJPT All right reserved
Research J. Pharm. and Tech 2017; 10(9):3095-3097.
DOI: 10.5958/0974-360X.2017.00549.2